Ab initio calculations for bromine adlayers on the Ag(100) and Au(100) surfaces: 

the c(2 x 2) structure 
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Ab initio total-energy density-functional methods with supercell models have been employed to 
calculate the c(2 x 2) structure of the Br-adsorbed Ag(100) and Au(100) surfaces. The atomic 
geometries of the surfaces and the preferred bonding sites of the bromine have been determined. 
The bonding character of bromine with the substrates has also been studied by analyzing the 
electronic density of states and the charge transfer. The calculations show that while the four- 
fold hollow-site configuration is more stable than the two-fold bridge-site topology on the Ag(100) 
surface, bromine prefers the bridge site on the Au(100) surface. The one-fold on-top configuration 
is the least stable configuration on both surfaces. It is also observed that the second layer of the Ag 
substrate undergoes a small buckling as a consequence of the adsorption of Br. Our results provide 
a theoretical explanation for the experimental observations that the adsorption of bromine on the 
Ag(IOO) and Au(IOO) surfaces results in different bonding configurations. 



PACS numbers: 68.43.Bc, 68.43.Fg, 68.47.-b, 82.45.-h 



I. INTRODUCTION 

Anion adsorption on metals can strongly modify sur- 
face morphology and electronic structure and chemical 
reactivity. It is therefore of great scientific and techno- 
logical importance. In particular, the halide-adsorbed 
noble-metal systems play a significant role in electro- 
chemistry. From the fundamental point of view, halide- 
adsorbed noble-metal surfaces are important model sys- 
tems for adsorption on metal surfaces with formation of 
ordered two-dimensional adsorbate structures. It is thus 
not surprising that the adsorption of halidcs on noble 
metals has been extensively investigated. 

The systems we selected to study are the Br- 
chemisorbed Ag(100) and Au(100) surfaces. The adsorp- 
tion of bromine on the Ag(100) and Au(100) surfaces 
both in vacuum-and in solution have been widely-studied 
by experimentaJTl3 and by classical simulations .EjO Ex- 
perimentally, bromine has been found to form different 
bonding structures on the Ag(100) and Au(100) surfaces. 
While bromine chemisorbed on the Ag(100) surface occu- 
pies the four- fold hollow site (hereafter referred to as H4), 
the most stable chemisorption structure on Au(100) is 
the configuration with bromine at the two-fold bridge site 
(hereafter referred to as B2). These different chemisorp- 
tion structureSrhawe been verified by various experimental 
measurements. liralailS However, theoretical studies have 
not yet reproduced these different adsorption behaviors. 

Kleinherbers et al. performed angle-resolved photoe- 
mission, low-energy electron diffraction (LEED), and X- 
ray photoemission measurements for the interaction of 
halides with Ag surfaces n They found that the adsorp- 
tion of CI, Br, and I on the Ag(100) surface in vacuum all 
resulted in the formation of a c(2 x 2) overlayer with the 
adsorbates in the H4 sites. Using in situ surface X-ray 



scattering, Ocko et al. studied the adsorption of bromide 
on an Ag(IOO) electrode. They observed a disordered 
phase at lower coverages and an-.ordered c(2 x 2) phase at 
a coverage of half a monolayer .□ The Br was determined 
to bond at the H4 site in both the c(2 x 2) and disor- 
dered phases. The atomic geometry of the Br/Ag(100)- 
c(2 x 2) surface was further investigated by Endo et al. 
with the-.m situ X-ray absorption fine structure (XAFS) 
method.ta The H4 site was confirmed to be the bonding 
site of bromine. The disordered phase of the Br/Ag(I00) 
surface at lower bromine coverages was also recently fur- 
ther investigated experimentally.tj In that study, it was 
suggested that while most of the bromide ions occupy 
the H4 sites, there are additional bromide ions adsorbed 
slightly off the H4 sites. 

The LEED data reported by Bertel et alEB have shown 
that the chemisorption of Br on the Au(100) surface in 
vacuum results in the rearrangement of the top-layer Au 
atoms of the original clean reconstructed Au(100)-(5x 20) 
surface and the formation of an unreconstructed (lxl) 
substrate structure. Several ordered structures of the 
Br adlayer, including c(2 x 2), (y/2 x 4V2)R45°, and 
c(4 x 2)i?45°, were obtained after bromine exposure on 
Au(100) surfaces, with the former two structures being 
metastable. It was concluded from the experimental data 
that Br adsorbed at the B2 site on the Au(100) surface in 
all the observed phases. This is in contrast to the case of 
the Br-adsorbed Ag(100) surface. Under electrochemical 
in situ conditions, surface X-ray scattering and scanning 
tunneling microscopy (STM) experiments showed that 
bromide adsorbed on the unreconstructed Au(100)-(1 x 1) 
surface forms a commensurate c(v / 2 x 2\/2)-R45 struc- 
ture and an incommensurate c(\/2 x_2n)JM5 (p < 2v2, 
depending on the applied potential) BElOO In this case, 
too, the bromide ions were determined to reside at the 
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J5 2 sites. 

In contrast to the considerable progress of the exper- 
imental measurements, theoretical studies employing ab 
initio methods to these systems are still at an early stage. 
Several groups have performed ab initio Hartree-Fock 
(HF) and density-functional-theory (DFT) calculations 
for the Br/Ag£l00) and Br/Au(100) interfaces using clus- 
ter models. U~E3 While these investigations have provided 
useful information about the interaction between Br and 
the surfaces, as we discuss below, many of the results are 
not yet sufficiently accurate. For example, the preferred 
bonding site of Br on the Au(100) surface was incor- 
rectly predicted by ab initio DFT cluster calculations, 
which showed that Br would prefer to bond atjJhe H4 
site on both the Ag(100) and Au(100) surfacesiia Simi- 
larly, ab initio HF studies with small clusters predicted 
the bridseiSite as the preferred adsorption site for Br on 
Ag(100).lia Given that these calculations with small clus- 
ters cannot reproduce such a fundamental property as 
the binding site, all other results (e.g., energy barriers) 
obtained from such calculations for both the Br/Ag(100) 
and Br/Au(100) systems are questionable. 

Here we present results of total-energy DFT calcula- 
tions in which we used supercell models for the Ag(100) 
and Au(100) surfaces. The detailed atomic structures 
and electronic properties of the chemisorbed surfaces 
and the preferred bonding site of the adsorbate have 
been determined. Our theoretical approach has repro- 
duced the different behavior of Br on the Ag(100) and 
Au(100) surfaces. The most stable adsorption sites for 
Br chemisorbed on the Ag(100) and Au(100) surfaces are 
determined by our calculations to be the H4 and the B2 
sites, respectively. Our results are in excellent agreement 
with the experimental data. The obtained results for the 
electronic properties also enable us to analyze the nature 
of the bonding between Br and the substrates and un- 
derstand the different adsorption behavior of Br on the 
Ag(100) and Au(100) surfaces. 

The remainder of this paper is organized as follows. 
In Sec. II we outline in detail the computational method 
and the supercell models that we used. In Sec. Ill we 
present and discuss the results for bulk, clean surfaces 
(Sec. Ill A), and adsorbed surfaces (Sec. Ill B). The ad- 
sorption geometries and atomic relaxations are discussed 
in Sec. Ill Bl, and the electronic properties and bonding 
character in Sec. Ill B2. We also give comparisons of our 
results with previous calculations and experimental data. 
Finally, in Sec. IV, we summarize the main results of our 
calculations. 



II. METHOD AND MODEL 

On the Ag(100) and Au(100) surfaces, there are three 
different symmetric adsorption sites, known as H4, B2, 
and T\ (on-top) sites. These three sites are shown in 
Fig. 1. We have studied a c(2 x 2) structure in which Br 




(a) (b) 




(c) 

FIG. 1. Schematics of an adatom at (a) the four- fold hol- 
low (Hi) site, (b) the two-fold bridge (B2) site, and (c) the 
on-top (Ti) site on the unreconstructed Ag(100) and Au(100) 
surfaces. 



forms an adlayer with a coverage of i monolayer on the 
surface for each of the three bonding configurations. 

The metal surface is modeled by repeated slabs with 
five, seven, and nine metal layers separated by a vac- 
uum region equivalent to five or seven metal layers. Each 
metal layer in the supercell contains two metal atoms. Br 
is adsorbed symmetrically on both sides of the slab. All 
the metal atoms were initially located at their bulk po- 
sitions, with the equilibrium lattice constant of the bulk 
determined by our calculations. 

The calculations were performed within density- 
functional theory, using the pseudopotential (PP) 
method and a plane- wave basis set. The results re- 
ported in this paper were obtained using-, the Vi- 
enna ab-initio simulation package (VASP).ElrE^I The 
exchange-correlation effects were treated with the gener- 
alized gradient-corrected exchange-cor:piat|ion function- 
als (GGA) given by Perdew and Wang£a£3 We adopted 
the scalar-relativistic Vanderbilt ulteasoft pseudopoten- 
tials supplied by Kresse and HafherE3'EZl. A plane-wave 
energy cutoff of 20 Ry and 56 special k points in the ir- 
reducible part of the two-dimensional Brillouin zone of 
the c(2 x 2) surface were used for calculating both the 
Br/Ag(100) and Br/Au(100) surfaces. Optimization of 
the atomic structure was performed for each supercell 
via a conjugate-gradient technique using the total energy 
and the Hellmann-Feynman forces on the atomsBfl All 
the structures were fully relaxed until the change in total 
energy was smaller than 1 meV between two ionic steps. 
The convergence of the total energies was checked with 
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different values of the plane-wave cutoff and different 
numbers of special k points. A series of test calculations 
with different slab thicknesses (from five to nine metal 
layers) and vacuum-gap widths (equivalent to five and 
seven metal layers) were also carried out to check con- 
vergence. The calculations on which this paper is based 
represent approximately 300 CPU hours on an IBM SP2 
computer. 

III. RESULTS AND DISCUSSION 
A. Bulk and clean surface 

We first present the calculated properties for bulk sil- 
ver and bulk gold, and the relaxed but unreconstructed 
clean Ag(100) and Au(100) surfaces. 

Calculations for bulk Ag and Au were conducted with 
408 special k points and cutoff energies ranging from 20 
Ry to 40 Ry. The total energy convergence with respect 
to the cutoff energy was shown to be within a few tenths 
of 1 meV. We obtained lattice constants of 4.17 A and 

4.18 A for bulk Ag and Au, about 2.0% and 2-5% larger 
than the corresponding experimental valuesEa at room 
temperature, respectively. Previous total-energy DFT 
calculations at the GGA level found la|ttice, constants be- 
tween 4.13 A and 4.19 A for bulk Ag,ErE3 and between 

4.19 A [Rcf. H and 4.20 A [Ref. |7| for bulk Au. Our 
results are in good agreement with these calculations. 

The properties of the clean Ag(100) and Au(100) sur- 
faces were calculated using supercells containing a 7-layer 
metal slab and a vacuum gap equivalent to 7 bulk-metal 
layers. A 1 x 1 surface cell was used with 66 special 
k points in the surface Brillouin zone. The kinetic en- 
ergy cutoff for the calculations was 20 Ry. All the layers 
except for the central one were relaxed. Surface recon- 
struction was not considered. 

The surface energies for the Ag(100) and Au(100) sur- 
faces obtained from our calculations are 0.43 eV/atom 
and 0.47 eV/atom, respectively. Both results are in 
good agreement with recent pscudopotential GGA cal- 
culations by Yu and Scheffler, who reported—the corre- 
sponding values of 0.48 and 0.45 eV/atom.Oo A re- 
cent calculation with linear-muffin-tin-orbital (LMTO)- 
GGA methods, however, obtained much larger values of 
0.65 and 0.90 eV/atom for the. unrelaxed Ag(100) and 
Au(100) surfaces, respectively!!!] The reason for the large 
discrepancy from the other GGA results is not cigar, 
PseudopotentiaUocal-density approximation (LDA) ,EjE£| 
LMT^-LDAj2jii2l and linearized augmented-plane-wave 
LDAcjCj calculations have provided values of 0.59-0.7 
eV/atom, and 0.69-0.72 eV/atom for the surface ener- 
gies of the Ag(100) and Au(100) surfaces, respectively. 
These DFT-LDA values for the surface energy are gen- 
erally larger than those calculated with the DFT-GGA 
calculations reported in this paper and Refs. |3l] and |3"q . 
This is consistent with the previous observations that 



LDA surface energies are normally larger than the cor- 
responding GGA values due to the different treatment 
of the exchange-correlation functional. The calculated 
values for the surface energies are thus seen to be quite 
sensitive to the computational method and the form of 
the exchange-correlation functional. 

Table I shows the results of the surface relaxation. 
While no significant structural relaxation is found for ei- 
ther surface, the Ag(100) surface shows a slightly more 
relaxed geometry than the Au(100) surface. Both sur- 
faces show an inward relaxation of the top layer and 
slight outward relaxation .of the second and third lay- 
ers. LEED measurements^ showed insignificant relax- 
ation of the Ag(100) surface with Adi 2 /d = ± 1.5% 
and Ad23/do — ± 1.5%, where Ae?i2 and A<i23 are the 
changes in spacing between the top and the second layer 
and between the second and the third layer, and do is 
the bulk interlayer distance. Our results are thus in good 
agreement with the experimental data, and basically con- 
sistent with other ab initio calculations, the results of 
which are also listed in Table I for comparison. 

B. Br-adsorbed Ag(lOO) and Au(lOO) surfaces 

1. Relaxations and energetics 

The results of our calculations for the Br/Ag(100)- 
c(2 x 2) and the Br/Au(100)-c(2 x 2) surfaces are shown 
in Tables II and III. If not otherwise indicated, the re- 
sults reported in this section (and in Tables II and III) 
were calculated with a cutoff energy of 20 Ry, 56 special 
k points, and supercells containing a 9-layer metal slab 
with a vacuum region equivalent to seven metal layers. 



TABLE I. Relaxation of the clean Ag(100) and Au(100) 
surfaces. Aalij is the change of the interlayer distance, and do 
is the corresponding distance in the bulk. 





Adi 2 /do (%) 


Ad 23 /d (%) 


Ad 34 /d (%) 


Ag(100) 
This work 


-1.8 


0.7 


0.2 


PP-GGA a 


-1.4 






PP-LDA a 


-2.2 


0.4 




PP-LDA b 


-1.3 


1.0 


0.8 


LMTO-LDA c 


-1.9 






Experiment* 1 


0±1.5 


0± 1.5 




Au(100) 
This work 


-1.3 


0.3 


0.2 


PP-LDA C 


-1.2 


0.4 




LMTO-LDA f 


-1.0 







a Ref. |T 
b Ref. M 
c Refs.J|-|l[ 
d Ref. 0a 
c Ref. Bsf 
f Refs.^g,|l|. 
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In-plane relaxations of the top-layer metal atoms were 
found to result in changes of the distance between Br 
and its nearest-neighbor metal atoms and of the total- 
energy difference between two different configurations 
within only 0.01 A and a few meV, respectively. The ef- 
fects of in-plane relaxations are thus negligible, and such 
relaxations were not considered in the calculations for the 
Br-adsorbed surfaces. 

Table II shows total-energy differences between differ- 
ent bonding configurations of both surfaces. Each struc- 
ture was optimized. We found that while the total energy 
of the H4 configuration is lower by 213 meV than the B 2 
configuration for the Br/Ag(100) surface, it is higher by 
58 meV for the Br/Au(100) surface. The T x configura- 
tion for both surfaces is found to be higher in total energy 
than both the corresponding H4 and B 2 configurations. 
Thus, we conclude that while Br adsorbed on the Ag(100) 
surface prefers the H4 site, it is adsorbed at the B 2 site 
on the Au(100) surface. This coadusioji is in agreement 
with experimental observations.EJcl'ErEj 

It is interesting to note that the magnitude of the total- 
energy difference between the H4 and B 2 structures for 
Br/Ag(100) is significantly larger than the correspond- 
ing value for the Br/Au(100) surface. This suggests that 
diffusion of Br on the Au(100) surface may occur much 
more easily than on the Ag(100) surface since the total- 
energy difference between the most stable configuration 
(the global minimum) and the less favorable configura- 
tion (probably a saddle point) is directly relevant to ad- 
sorbate diffusion. 

Previous theoretical studies employing cluster mod- 
els also determined the preferred bonding sites of Br 
on the Ag(100) and Au(100) surfaces, as mentioned in 
Sec. I. Ah initio HF calculations showed that Br would 
prefer to bond at the B 2 site on the Ag(100) surface 
(by 370 meV/adatom over the H4 site, and by 570 
meV/adatom over the T\ site).Ej This is inconsistent with 
both our—DFT-supercell calculations and the experimen- 
tal dataEcl DFT cluster calculations predicted that the 
binding energy of Br at the H4 site on both the Ag(100) 
and Au(100) surfaces was larger than at the B 2 and Tf 
sites by 120 meV for Ag(100) [89 meV for Au(100)] and 



TABLE II. Total energy differences (in eV per unit cell) 
between different configurations of the Br/Ag(100)-c(2 x 2) 
and Br/Au(100)-c(2 x 2) surfaces, obtained from calculations 
with supercells containing a 9-layer slab and a 7-layer vacuum 
region, a cutoff energy of 20 Ry, and 56 special k-points. 



Eh 4 


-E B2 (Br/Ag(100)) 


-0.213 


Er 4 


- E Tl (Br/Ag(100)) 


-0.557 


Eb 2 


- E Tt (Br/Ag(100)) 


-0.344 


Eh 4 


-E B2 (Br/Au(100)) 


+0.058 


Eh 4 


-E Tl (Br/Au(100)) 


-0.244 


Eb 2 


-E Tt (Br/Au(100)) 


-0.302 



202 meV for Ag(100) [202 prefers to bond at the H 4 site 
on the Au(100) surface, is-.in disagreement with exper- 
imental measurements, Bii3 as well as with our results. 
We believe that the main problem is that these previ- 
ous calculations were limited to small clusters, contain- 
ing only up to 13 metal atoms. It is well known that a 
small metal cluster has a very different electronic struc- 
ture than an extended metal surface, yielding very sig- 
nificant differences. in.adsorbate binding energies and re- 
action pathwaysoTJ Large clusters or extended surface 
models (e.g., supercell models) are therefore needed to 
simulate real metal surfaces accurately. 

The structural parameters of the optimized geometries 
for the H4, B 2 , and T\ configurations of the Br/Ag(100) 
and Br/Au(100) surfaces are presented in Table Iff. The 
vertical distances (d z ) between the Br centers and the 
plane of the centers of the top-layer atoms were calcu- 
lated to be 1.91 A on the Ag(100) surface and 2.01 A 
on the Au(100) surface for the H4 structure, 2.16 A on 
Ag(100) and 2.18 A on Au(100) for the B 2 configuration, 
and 2.48 A on Ag(100) and 2.46 A on Au(100) for the 
T\ structure. While the values of d z for the B 2 and T\ 
configurations of the Br/Ag(100) surface are very close 
to the corresponding values for the Br/Au(100) surface, 
the distance between Br and the surface in the H4 config- 
uration is observed to be significantly longer (by 0.1 A) 
on Br/Au(100) than on Br/Ag(100). Accurate in-situ 
XAFS measurements for Br/Ag(100) in NaBr solution 
by Endo et al. showed that the bond length between Br 
and its four nearest-neighbor Ag atoms in the H4 con- 
figuration is 2.82 ± 0.05 A, and the .distance between the 
Br and the surface is f .94 ± 0.07 A.l Our results (2.82 A 
and 1.91 A, respectively) are thus in excellent agreement 
with the experimental data, provided that the solution 
has only a minor influence on the bond lengths between 
the adsorbate and the surface. 

The bond lengths between Br and Ag and Au clusters 
of varying size have been obtained with both HF and 
DFT calculations. Illas et al., using the HF method with 
a cluster of 5 Ag atoms simulating the H4 configuration 
of the Ag(100) surface, obtained a value of 3.43 A for 
the length of the Br-Ag bond.E3 Paccioni, also using the 
HF method with slightly larger clusters, found that the 
bond lengths between a Br ion and the surface Ag atom 
were 3.24 A, 2.97 A, and 2.94 A in clusters of Br~-Ag 13 
(modeling the H4 structure), Br~-Ag 8 (simulating the 
B 2 geometry), and Br~-Agi3 (representing the Xi con- 
figuration) , respectively, fgnaczak and Gomes performed 
DFT calculations with clusters containing a Br ion and 
12 metal atoms and determined the bond lengths to be 
3.2 A, 3.0 A, and 2.9 A for the H4, B 2 , and T x configu- 
rations of Br~-Agi2 and Br~-Aui2 clusters, respectively. 
All of these values are much larger that those obtained 
from our supercell calculations and the XAFS measure- 
ments, suggesting that small clusters do not represent 
the metal surfaces properly. 

The B 2 configuration of the Br/Ag(100) surface shows 
a very similar relaxation of the surface metal layers as 
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that of the same configuration for the Br/Au(100) sur- 
face. Both undergo an inward relaxation of the top layer 
and slight outward relaxations of the second and third 
layers. Similar relaxed structures are also found for the 
clean Ag(100) and Au(100) surfaces (see Table I). 

Our calculations show that the second metal layer un- 
dergoes a small buckling with the adsorption of Br in the 
Hi configurations. The atoms in the second layer that 
are immediately below the H4 sites are observed to shift 
slightly up towards the surface, while the other atoms in 
that layer shift up by only on the order of 0.001 A and 
hence essentially keep their bulk positions. The spacing 
between these two sub-layers is found to be 0.02 A and 
0.04 A for the Br/Ag(100) and Br/Au(100) surfaces. The 
distance between Br and the second-layer metal atom 
just below it is still far larger than the bond length be- 
tween Br and its nearest-neighbor metal atoms in the 
top layer. Thus a pseudo-five-fold coordination, which 
has been obserited-in c(2 x 2) overlayer structures on bcc 
metal surfaces, lIe3 does not exist for the Br/Ag(100) and 
Br/Au(100) surfaces. This buckling may give rise to an 
effective Br-Br interaction, mediated through the surface 
strain field. The top metal layer in the H4 configurations 
still shows a slight inward relaxation, similar to the cases 
of the clean surfaces and the B2 configurations. 

The top metal layer of the T± configuration of the 
Br/Au(100) surface also shows a small buckling. The 
Au atoms in the top layer that are bonded to Br are 
observed to undergo a larger inward relaxation than the 
other half of the Au atoms in the top metal layer. The 
corresponding buckling is, however, very large for the T\ 
configuration of the Br/Ag(100) surface. The distance 
between the two sublayers formed from the top Ag layer 
is 0.75 A, indicating a zigzag surface reconstruction. 

Finally, in Table IV we show results of convergence 
checks for the total-energy differences. Such checks are 
particularly important for the Br/Au(100) surface due to 
the small value of the total-energy difference between the 
Hi and B2 configurations. Calculations with a higher 
cutoff energy (30 Ry) obtained total-energy differences 
within 1 meV of those from calculations with a cutoff 
energy of 20 Ry. Increasing the number of special k 
points from 36 to 56, increasing the slab thickness from 
7 to 9 metal layers, and increasing the vacuum region 
in the supercell from 5 to 7 layers, all changed the re- 
sults by only a few meV. Supercells with 5 metal layers 
are seen to cause errors in the total-energy differences of 
~30 meV for Br/Ag(100) and -10 meV for Br/Au(100). 
The use of 20 special k points also causes an error of 
—10 meV. Therefore, it is necessary to employ supercells 
with at least 7 metal layers and 36 special k points for 
obtaining the total-energy differences with errors smaller 
than 10 meV. The distances between Br and its nearest- 
neighbor metal atoms were also checked. We found that 
the changes of these distances were smaller than 0.01 A 
over the ranges of cut-off energies between 20 and 30 Ry, 
numbers of k points between 20 and 56, and numbers of 
metal layers between 5 and 9 in the supercells, indicating 



that the bond lengths are not very sensitive to the choice 
of computational parameters. 



2. Electronic properties and bonding character 

In order to better understand the differences between 
the bonding of bromine on the Ag(100) and Au(100) sur- 
faces, we calculated the total electronic density of states 
(DOS), the DOS projected onto individual atoms and 
specific atomic states, and the charge transfer between 
bromine and the substrate. 

Figure 2 shows the total DOS for the Br-adsorbed 
Ag(100) and Au(100) surfaces. For comparison, the total 
DOS for the clean relaxed Ag(100) and Au(100) surfaces 
are also shown. The peaks of the DOS curves for the 
clean surfaces represent the main features of the s and 
d states of the substrates and remain essentially at the 
same positions when bromine atoms are adsorbed. New 
states are, however, found to be located at between —15 
eV and —13 eV relative to the Fermi level in the DOS 
curves of the Br-adsorbed Ag(100) and Au(100) surfaces. 
These states arc predominantly the bromine 3s states 
with small contributions from the s and d states of the 
substrate, as seen from the curves for the DOS projected 
onto the specific atomic states of the adsorbate and the 
substrate (shown in Fig. 3). Significant changes of the 
total DOS in the higher energy range close to the Fermi 
level (above -2.5 eV and -1.5 eV for the Br/Ag(100) and 
the Br/Au(100) surfaces, respectively) are also observed 
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FIG. 2. Total density of states for the clean relaxed and 
the Br-adsorbed Ag(100) and Au(100) surfaces. The Fermi 
level is at eV. In this figure, as well as in Figs. 3 and 
4, the curves are obtained from calculations with supercells 
containing a 7-layer slab and a 7-layer vacuum region, a cutoff 
energy of 20 Ry, and 36 special k-points. 
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FIG. 3. The density of states projected onto the Br 3s, 
Br 3p, Ag 5s, Ag 4d, Au 6s, and Au 5d states for the H4, 
B2, and Xi configurations of the Br-adsorbed Ag(100) and 
Au(100) surfaces. 



when the DOS curves for the clean surfaces are compared 
with those for the H4, B2, and T\ configurations of the 
Br-adsorbed surfaces (see Fig. 2). The electronic states 
in the higher energy range are composed mostly of the 
bromine 3p states and the Ag Ad (or Au 5d) states, with 
some contributions coming from the Ag 5s (or Au 6s) 
states (see Fig. 3). 

In Fig. 4, we show the DOS projected onto bromine for 
the systems before and after adsorption. The results for 
the systems before adsorption were calculated by employ- 
ing the supercell of a 7-layer slab and a 7-layer vacuum 
region. A bromine layer with a c(2 x 2) periodicity was 
kept fixed in the middle of the vacuum region (located 
at ~ 7 A above the surface) so that there was essen- 
tially no interaction between Br and the substrate. The 
peaks located in the lower and higher energy ranges in 
the projected DOS before adsorption are due to the Br 



3s and 3p states, respectively. The slight broadening of 
the 3p states reflects the weak 3p-3p interaction between 
neighboring bromine atoms. When bromine is adsorbed, 
both the 3s and 3p states shift down in energy due to the 
bonding between bromine and the substrate. A broad- 
ening of the Br 3p states is also observed and can be 
attributed to the hybridization of the bromine 3p states 
with the s- and d-bands of the substrates (see Fig. 3). 
The Br 3s states are also seen to mix slightly with the 
s and d states of the substrate (see also Fig. 3). The 
hybridization of the Br 3s and 3p states with the elec- 
tronic states of the substrate suggests covalent bonding 
between bromine and the Ag(IOO) and Au(IOO) surfaces. 

The bonding of Br with the Ag(100) and Au(100) sur- 
faces is also found to be associated with a charge transfer 
from the substrate to Br. To obtain a rough estimate 
of the charge transfer, we calculated the change of the 
charge for a bromine atom upon adsorption by integrat- 
ing the difference of the corresponding charge densities, 
over a sphere with a radius of 1.28 A around the atom.L^I 
We found that 0.15 and 0.14 electrons were transfered 
from the Ag(100) and Au(100) surfaces, respectively, to 
the bromine atom. The amount of the charge transfer 
was found to be basically the same for the H4, B2, and T\ 
configurations. These results are consistent with the data 
for the DOS projected onto the Br atom. By integrating 
the 3p contributions up to the Fermi level, we observe 
that more 3p states are occupied in the Br-adsorbed sur- 
faces than in the systems before adsorption (see Fig. 4). 
A recent periodic GGA calculation with a local basis 
set for the adsorption of chlorine on the Ag(lll) sur- 
face also found that a slight charge (~ 0.2 electrons) was 
transferred from the Ag(lll) substrate to the chlorine 
atom.Lj Experimental measurements of the electrosorp- 
tion valency of Br adsorbed on p Ag j jl QQj report values 
of approximately —0.70 to — 0.75,Erli3tJ'Fl3 corresponding 
to a residual charge of 0.25 to 0.30 electrons on the ad- 
sorbed Br. These values are considerably larger than our 
calculated charge of 0.15 electrons. The discrepancy may 
be due to the fact that our calculations were performed 
for systems in vacuum. In an electrochemical environ- 
ment, the net charge associated with the adsorbate might 
be very different from that in vacuum, due to solvation. 
However, the discrepancy might also be attributed to in- 
accuracies in the theoretical and experimental methods 
used to estimate the charge. 

The difference in the bonding strength of the bromine 
with the substrate between the different configurations 
directly affects their relative stability. Based on our DOS 
data, we provide a qualitative explanation of the differ- 
ence in bonding strength between the H4 and B2 config- 
urations. The Br 3s and 3p states in the Br/Ag(100)-i/4 
configuration are significantly lower in energy than in the 
Br/Ag(100)-i?2 configuration (see Fig. 4). In addition, 
the intensity of the lower part of the 3p states (below ap- 
proximately —3 eV) is larger for Br/Ag(100)-fZ4 than for 
Br/Ag(f 00)-B 2 - Both facts suggest a stronger bonding 
for the H4 configuration on Ag(100). This is expected 



G 



since there are more direct bonding neighbors for the 
bromine atom at the H4 site. On the other hand, the Br 
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FIG. 4. The density of states projected onto Br for the H4, 
B2, and Ti configurations of the Br-adsorbed Ag(IOO) and 
Au(IOO) surfaces (a) over a larger energy range (—16 eV to 
2 eV) in which both the 3s and 3p states are shown, and (b) 
over a smaller energy range (—8 eV to 2 eV) where only the 3p 
states are presented. Also shown in (a) is the DOS projected 
onto Br without adsorption. 



3s states in the Br/Au(100)-_ff4 configuration are only 
slightly lower in energy than in the Br/Au(100)-i?2 con- 
figuration. While the Br 3p states extend over almost 
the same range in energy for the Br/Au(100)-£/4 and - 
B2 configurations, they have slightly larger intensity in 
the lower part (below ~ —2 eV) and smaller intensity in 
the higher part (above ~ -2 eV) for the Br/Au(100)-i?4 
configuration than for the Br/Au(100)--B2 configuration. 
While the H4 configuration thus has a stronger covalent 
bonding for both the Br-adsorbed Ag(100) and Au(100) 
surfaces, the difference in bonding strength between the 
H4 and the corresponding B2 configurations is smaller for 
the Au(100) surface than for the Ag(100) surface. This is 
probably due to the fact that the Au 6s and 5d electrons 
are more delocalized than the Ag 5s and Ad electrons, 
and the bonding strength is expected to be less sensitive 
to the bonding sites for substrates with more delocalized 
electrons. In addition to the stronger covalent bonding, 
the Coulomb attraction resulting from the charge transfer 
in the H4 configuration is also stronger than the corre- 
sponding B 2 configuration for both the Br/Ag(100) and 
Br/Au(100) systems, due to the shorter distance between 
bromine and the surface in the H4 configuration (see Ta- 
ble III). 

It is clear that there is a delicate competition between 
the attractive and repulsive interactions in each config- 
uration. In particular, the core-core repulsion between 
bromine and the substrate, which is irrelevant to the elec- 
tronic DOS but makes a contribution to the total energy 
of the system, in the H4 configuration is stronger than 
in the B2 configuration. The coxe=core repulsive energy 
is calculated as an Ewald sumE3i£3 (see the third col- 
umn of Table V) . The total energy, as determined in our 
DFT calculations, contains as separate parts the elec- 
tronic and the core-core Coulomb contributions. Differ- 
ences of the electronic and the core-core contributions to 
the total energy between the H4 and the B2 configura- 
tion are presented in Table V. We note that considera- 
tion of the electronic contributions alone does not prop- 
erly address the opposite order of the total-energy dif- 
ference (ie, the binding-energy difference) between the 
H4 and B2 configurations for the the Br/Ag(100) and 
Br/Au(100) surfaces. The core-core interactions need to 
be included. For both the Br/Ag(100) and Br/Au(100) 
surfaces, the electronic contribution favors the H4 con- 
figuration, while the core-core contribution favors the B2 
configuration. In the Ag(100) case, the core-core energy, 
which is higher for H4 than for B2 , is more than compen- 
sated by the lower electronic energy for the H4 configu- 
ration, resulting in H4 being the preferred bonding site 
for Br on Ag(100). For the Br/Au(100) system, how- 
ever, the lower electronic energy for the H4 configuration 
only partially compensates the higher core-core energy 
for Br/Au(100)-if4. As a result, the B2 configuration is 
lower in total energy than the H4 configuration for the 
Br-adsorbed Au(100) surface. The small magnitudes of 
the total-energy differences, compared to the individual 
electronic and core-core contributions, strongly cmpha- 
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size the need for very accurate energy calculations and 
careful convergence checks. 



IV. CONCLUSIONS 

The theoretical approach of supcrccll models combined 
with first-principles total-energy DFT pscudopotcntial 
methods has reproduced experimental measurements of 
preferred adsorption sites for Br-chemisorbed Ag(100) 
and Au(100) surfaces. 

We have shown that while the hollow-site configu- 
ration is more stable on the Br/Ag(100) surface (by 
210 meV/adatom over the bridge-site structure), the 
bridge-site configuration is more stable than the corre- 
sponding hollow-site structure by 60 meV/adatom on 
the Br/Au(100) surface. The calculations also predict 
that the one-fold on-top configuration is the least sta- 
ble structure on both surfaces (560 meV and 300 mcV 
higher than the corresponding most stable structure for 
the Br/Ag(100) and Br/Au(100) surfaces, respectively). 
Other aspects of the geometries of the Br/Ag(100) and 
Br/Au(100) systems have also been determined and are 
shown to be in excellent agreement with the available 
experimental data. 

The bond between Br and the substrate is found to be 
covalent with a slight polarization due to a small charge 
transfer from the substrate to the bromine. The chemi- 
cal bonding between Br and the substrate is shown to be 
stronger in the H4 configuration than in the B 2 configu- 
ration. Compared with the Br/Ag(100) surface, however, 
the Br/Au(100) surface exhibits a reduced difference in 
the bonding strength between the H4 and B2 configu- 
rations. The core-core Coulomb interaction is found to 
be higher for the H4 configuration than for the B2 con- 
figuration. The detailed balance between the electronic 
and the core-core contributions to the total energy de- 
termines H4 and £>2 as the the preferred bonding site on 
the Ag(100) and Au(100) surfaces, respectively. 

Our work demonstrates that the use of extended sur- 
face models and careful convergence checks are critical 
for obtaining reliable information on the Br/Ag(100) and 
Br/Au(100) systems from ab initio calculations. 
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TABLE III. Relaxation of the H 4 , B 2 , and Ti configurations of the Br/Ag(100)-c(2 x 2) and Br/Au(100)-c(2 x 2) surfaces. 
Also shown are the vertical distance (d z , in units of A) between Br and the plane of the centers of the top-layer atoms on the 
surface, and the distance (d, in units of A) between Br and its nearest- neighbor metal atom(s). The change in spacing between 
layers i and j is denoted by Adij. When sub-layers are present, the changes are denoted by Adij and Ad^, with Adij being 
the larger in magnitude, do is the same as in Table I. The computational details are the same as in Table II. 





d z 


d 


Ad 12 /d a (%) 


Ad' 12 /d (%) 


Ad 23 /d (%) 


Ad' 23 /d (%) 


Ad 34 /d (%) 


Br/Ag(100)-# 4 


1.91 


2.82 


-0.6 


0.1 


0.8 


0.1 


-0.0 


Br/Ag(100)-B 2 


2.16 


2.61 


-0.9 




0.5 




0.2 


Br/Ag(100)-Ti 


2.48 


2.48 


10.4 


-7.7 


0.2 




-0.1 


Br/Au(100)-# 4 


2.01 


2.89 


-1.6 


0.2 


1.4 


-0.3 


0.4 


Br/Au(100)-B 2 


2.18 


2.62 


-0.9 




0.4 




0.3 


Br/Au(100)-Ti 


2.46 


2.46 


-1.2 


-0.5 


0.5 




0.4 



TABLE IV. Convergence checks for the total energy differences (in units of eV per unit cell) 
between different configurations with respect to the cutoff energy (E cut , in units of Ry), the number 
of metal layers (N m ) in the supercell, the thickness of the vacuum region (N v , in units of number 
of bulk metal layers), and the number of k points in the surface Brillouin zone (iV fe ). 



Ecut 




N v 


N k 


Eh 4 — Eb 2 
Br/Ag(100) 


Eu A — Et x 
Br/Ag(100) 


Eh 4 — Eb 2 
Br/Au(100) 


Eb 2 — Et x 
Br/Au(100) 


20 


5 


7 


20 


-0.191 


-0.526 


+0.061 


-0.282 


20 


7 


7 


20 


-0.222 


-0.556 


+0.048 


-0.285 


20 


7 


7 


36 


-0.211 


-0.552 


+0.056 


-0.291 


30 


7 


7 


36 


-0.212 


-0.552 


+0.057 


-0.291 


20 


7 


5 


36 


-0.210 


-0.550 


+0.057 


-0.289 


20 


7 


7 


56 


-0.210 


-0.551 


+0.059 


-0.294 


20 


9 


7 


56 


-0.213 


-0.557 


+0.058 


-0.302 



TABLE V. Differences (in eV per unit cell) of two energy 
contributions to the total energy between the H4 and B 2 
configurations. E e and E cc are the electronic and core-core 
Coulomb contributions, respectively. Also shown is the to- 
tal-energy difference (E^* — Eg*). The results are obtained 
from calculations with supercells containing a 7-layer slab and 
a 7-layer vacuum region, a cutoff energy of 20 Ry, and 36 spe- 
cial k-points. 





rpe rpe 

E Hl - E B2 


^H 4 - ^B 2 


rptot rptot 
- &B 2 


Br/Ag(100) 


-429.230 


+429.018 


-0.212 


Br/Au(100) 


-371.156 


+371.213 


+0.057 
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